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Abstract 

The  Gamble  II  plasma-filled  rod  pinch  (PFRP) 
produces  an  order-of-magnitude  higher  x-ray  dose  and 
smaller  FWHM  source  size  than  achieved  with  a  vacuum 
rod  pinch.  Intense  electron  deposition  at  the  tip  produces 
a  high-energy-density  expanding-tungsten  plasma,  the 
bremsstrahlung  from  which  creates  a  low-intensity  halo 
around  the  source  core  that  increases  the  effective  source 
size  for  some  applications.  The  plasma  motion  has  been 
measured  with  a  holographic  interferometer.  Using  the 
distribution  of  electron  deposition  derived  from  line- 
spread  measurements,  the  observations  are  compared  with 
zero-and  one-dimensional  self-similar  MHD  modeling  of 
the  plasma  motion.  Results  demonstrate  that  magnetic- 
field  and  ohmic-heating  effects  during  electron  deposition 
are  important  to  understanding  PFRP  dynamics. 

I.  INTRODUCTION 

The  rod-pinch  diode  is  an  important  new  radiography 
source  in  the  1-MV  and  above  x-ray  regime. [1-5]  A 
standard  vacuum  rod-pinch  on  the  NRL  Gamble  II 
generator  [1],  consisting  of  a  mm-diam  tapered  tungsten 
anode  rod  protruding  through,  and  a  few  cm  beyond  an 
annular  carbon  cathode,  is  characterized  by  a  30-  to  60- 
Ohm  impedance  and  an  x-ray  source  size  comparable  to 
the  rod  diameter.  The  plasma-filled  rod  pinch  (PFRP) 
extends  megavolt  rod-pinch  radiography  to  low-imped- 
ance/high-current  generators  by  surrounding  the  rod 
downstream  of  the  cathode  with  a  1016-cm"3  plasma.  [6-7] 
The  plasma  initially  shorts  the  diode,  and  then  opens  to  a 
few-Ohm  impedance,  concentrating  a  0.5-MA,  MeV 
electron  beam  near  the  rod  tip.  The  high  current  density 
on  the  tapered  tip  produces  an  order-of-magnitude  higher 
x-ray  dose  (20  rad  at  1  m)  and  smaller  FWHM  source  size 
than  achieved  with  a  vacuum  rod  pinch  on  Gamble  II. 

The  intense  electron  deposition  in  the  PFRP  produces  a 
high-energy-density  expanding-tungsten  plasma,  the 
bremsstrahlung  from  which  creates  low-intensity  wings 
around  the  source  core  that  increases  the  effective  source 
size  for  some  applications.  This  two-dimensional,  time- 
dependent  tungsten  plasma  motion  has  been  studied  [7] 
with  a  holographic  interferometer  [8].  Using  an  axial 
distribution  of  electron  deposition  derived  from  side¬ 
viewing  (90°  to  the  rod)  line-spread  measurements  [5,9], 
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these  observations  are  compared  with  zero-  and  one¬ 
dimensional  self-similar  MHD  modeling  of  the  plasma 
motion.  Simulated  0°  (forward  on  the  rod  axis)  line 
spreads  derived  from  the  zero-dimensional  MHD  model 
agree  well  with  measurements.  Agreement  between  the 
one-dimensional  model  and  Schlieren  measurements  [7] 
demonstrates  that  magnetic-field  and  ohmic-heating 
effects  during  electron  deposition  are  important  to  under¬ 
standing  the  dynamics  of  the  plasma-filled  rod  pinch. 

II.  DISCUSSION  AND  CONCLUSIONS 

The  experimental  configuration  has  been  described 
elsewhere.  [6,7]  Here,  electrical,  Schlieren,  interfero¬ 
metry,  on-axis  and  side-viewing  rolled-edge  measure¬ 
ments  for  PFRPs  with  1 -mm-diam  tapered-tungsten  rods 
are  employed  in  the  analysis.  Figure  1  shows  the  diode 
current,  diode  voltage,  and  bremsstrahlung  x-ray  traces 
for  a  typical  shot.  Time  is  measured  from  the  start  of  the 
x-ray  pulse,  when  voltage  first  appears  at  the  rod  tip. 
Negative  times  correspond  to  the  run-down  phase  during 
which  the  injected  plasma  is  displaced  by  JxB  forces, 
sweeping  open  a  vacuum  gap  along  the  rod.  [6,7]  Voltage 
at  negative  times  is  then  associated  with  d(LI)/dt ,  where 
L  is  the  inductance  associated  with  the  evolving  vacuum 
gap.  Using  the  variations  of  V  and  I  shown  in  Fig.  1,  L  is 
calculated  to  be  close  to  linearly  rising  from  the  start  of 
the  current  pulse,  reaching  22  nH  at  t  =  0.  This  variation 
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corresponds  to  the  opening  of  a  cylindrical  volume  around 
the  rod  with  a  1 .2-cm  radius  and  length  increasing  from  0 
to  3  cm  during  the  run-down.  For  t  >  0,  the  voltage  at  the 
tip  Vtip  is  then  given  by  V  -  2.2x1  O'8 dl/dt,  so  that  the  elec¬ 
tron-beam  power  heating  the  full  3 -cm  length  of  the  rod  is 
taken  to  be  IVtiP  in  the  analysis. 

Figure  2  compares  Schlieren  images  of  vacuum  (50-kA) 
and  plasma- filled  (500-kA)  rod  pinches  60  to  70  ns  after 
the  start  of  x-radiation.  The  image  scales  are  identical  and 
the  rod  tips  are  aligned  (initial  rod  image  shown  for  the 
PFRP).  For  the  PFRP,  high-energy-density  tungsten 
plasma  expands  from  the  tip  during  bremsstrahlung  emis¬ 
sion,  creating  a  low-intensity  halo  around  the  x-ray-source 
core  that  is  manifested  as  extended  wings  in  the  on-axis 
line  spread  (Fig.  3).  Also  shown  in  Fig.  2  is  the  PFRP  90° 
line  spread  demonstrating  that  time-integrated  e-beam 
deposition  along  the  rod  correlates  with  the  observed 
plasma  expansion. 


Figure  2.  Schlieren  images  for  vacuum  and  plasma-filled 
RPs  60-70  ns  after  start  of  x-ray  emission. 


Figure  3  shows  typical  0°  edge-  and  line-spread  func¬ 
tions  for  the  PFRP  that  have  two  distinct  spot-size  scales. 
The  0.4-mm  FWHM  is  consistent  with  the  90°  line-spread 
distribution  of  Fig.  2  showing  concentration  of  electron 
deposition  on  the  tapered  tip.  The  1.7-mm  LLNL  defini¬ 
tion  source  size  [1]  reflects  the  extended  line-spread 
wings  of  Fig.  3,  which,  in  turn,  reflect  the  observed  rod- 
plasma  expansion  of  Fig.  2. 

The  zero-dimensional  MHD  analysis  seeks  to  explain 
the  observed  line-spread  wings  by  assuming  a  cylindrical 
self-similar  expansion  of  plasma  from  the  tapered  tip.  A 
Gaussian  tungsten-density  distribution  n(r,t)  of  width  R(t) 
is  assumed. 

n(r,t)=  N  -exp(-r2/R2)  (1) 

7iR2(t) 

Here,  N(l/cm)  is  the  fixed  line  density  given  by 

N  =  ikrRlp°  ’  (2) 

mw  is  the  mass  of  a  tungsten  ion,  p0  is  the  solid  rod  mass 
density,  Rb  =  L/20  is  the  radius  of  a  cone  of  rod-tip  length 
L  assumed  to  be  heated  by  the  e-beam.  Based  on  Fig.  2, 
L  «  3-4  mm  is  reasonable.  The  continuity  equation  is  then 


satisfied  when  the  radial  fluid  velocity  takes  the  form 


v(r,t) 


r  dR 
R(t)  dt 


(3) 


Substituting  these  forms  into  the  fluid  equation  of  motion, 

d2R  2(1 +  Z  )kT 

— —  =  ~ - .  (4) 

dt  mwR 


where  T(t)  is  the  temperature  and  Z(T)  is  the  ionization 
level.  Temperature  is  determined  from  radially-integrated 
energy  balance  per  unit  length  given  by 


Ph  -scrT  '2tiR 


\dt_rnwN_(dR 
r  2  \dt 


\2 


(5) 


In  Eq.  (5),  Eint  is  the  internal  energy,  the  last  term  repre¬ 
sents  kinetic  energy,  Ph  is  the  beam-heating  power  and 
black-body-like  radiation  Pbb  from  the  characteristic 
radius  R  is  assumed.  The  quantities  T(Eint)  and  Z(T)  are 
determined  from  SESAME  [10]  equation-of-state  data  for 
tungsten.  Measured  x-ray  dose  from  the  rod  tip  [6]  indi¬ 
cates  that  Ph  =  IV tip/ 3 L  is  appropriate;  a  piece-wise  linear 
fit  to  it  is  used  to  solve  Eqs.  (4)  and  (5).  Note  that  L  and 
emissivity  8  are  free  parameters. 

Figure  4  shows  solutions  of  the  above  equations  for 
8  —  0.3  and  two  values  of  L  =  3  and  4  mm.  For  each  trace, 
the  smaller  value  of  L  yields  the  upper  curve.  The  single 
Ph  curve  is  for  L  =  4  mm;  that  for  3  mm  is  4/3  larger.  At 
peak  temperature,  Z  in  the  14  to  16  range  is  calculated. 
Of  greatest  interest  here  is  R(t)  as  it  determines  the  wings 
of  the  0°  line  spread.  This  radius  is  found  to  depend 
weakly  on  s  in  the  range  0.03  to  0.3.  At  each  instant,  the 
bremsstrahlung  radiation  Pbr(r,t)  is  assumed  to  scale  like 
IVtip3/2  in  time  with  a  radial  dependence  proportional  to 
n(r,t)  of  Eq.  (1).  Integrating  Pbr  in  time  over  the  power 
pulse  then  provides  the  0°  point  spread  function  PS(r) 
used  to  calculate  the  line  spread  function  LS(y)  in  Eq.  (6). 
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Figure  4.  Zero-dimensional  tungsten-plasma  evolution 
for  L  =  3  mm  (upper)  and  L  =  4  mm  (lower). 

TC/  ,  f 00  PS ( r )rdr 

LS(y,-ljr 7  <*> 

Figure  5  compares  calculated  line  spreads  for  L  =  3  and 
4  mm  with  that  measured  in  Fig.  3.  Results  for  s  =  0.03 
are  indistinguishable  from  those  of  Fig.  5  for  e  =  0.3. 
These  good  fits  to  the  measurements  using  L  values  con¬ 
sistent  with  tip  deposition  from  the  90°  line  spread  (Fig.  2) 
confirm  that  PFRP  tip  expansion  is  the  source  of  the  ex¬ 
tended  wings  observed  on  the  0°  line  spread. 


position  (mm) 


Figure  5.  Measured  (points)and  predicted  (solid)  0°  line 
spreads  for  L  =  3  mm  (upper  curve)  and  L  =  4  mm  (lower 
curve). 


The  ability  of  the  zero-dimensional  model  to  predict  tip 
expansion  is  demonstrated  by  the  fits  of  Fig.  5,  and  con¬ 
sistency  between  the  values  of  R(t)  at  70  ns  in  Fig.  4  with 
the  Fig.  2  Schlieren  boundary  at  the  tip.  It  is  desired  to 
extend  the  analysis  to  a  one-dimensional  axial  self-similar 
model  in  order  to  better  understand  rod-plasma  expansion 
away  from  the  tip.  Equation  (1)  is  generalized  to  include 
the  axial  variation  of  line  density  associated  with  the  local 
solid-density  rod  radius. 

N( z)  =  nR2rod  ( z ) Pq  /  mw  (7) 


Local  rod  heating  Phz(z,t)  follows  the  90°  line  spread  of 
Fig.  2  in  space  and  IVtip  in  time.  The  local  beam  current 
density  Jbr(z,t)  entering  the  rod  plasma  follows  this  line 
spread  in  space  and  I  in  time.  Conservation  of  charge 
requires  that  the  axial  return  current  in  the  rod  Iz(z,t)  is 


Iz  =2tt\  ZgJbr(z'  ,t)rdz’  ,  (8) 


which  is  therefore  proportional  to  the  90°  edge  spread. 

With  the  addition  of  magnetic  pressure,  Eq.  (4)  is 
modified  to 


d2R(z) 

dt2 


2 

mwR 


(1  +  Z  )kT  - 


200N 


(9) 


in  cgs  units  except  for  Iz  in  A.  Energy-balance  now  takes 
the  form 


+  7.4x10 


6  ry (T)lj 
2jtR2 


1 1  dR 
100R  dt 


\ 

■  dt 


mwN  (  dR  V 

2  [~dt  J 


(10) 


where  rj(Ohm-cm)  is  the  Spitzer  resistivity.  The  ohmic 
heating  and  work  on  the  plasma  in  Eq.  (10)  derive  from 
the  radial  dependence  of  return-current  density  that  pre¬ 
serves  self-similarity  in  Eq.  (9),  i.e.,  jzxBe~  n(r,t). 


j2(r,z,t) 


KM  d 
2ttR2  udu 


[l-(l  +  u2)exp(-u2f2 


(11) 


In  Eq.  (11),  u  =  r/R(z,t).  The  self-similar  axial-current 
density  resembles  a  Gaussian  with  characteristic  radius 
1.2R.  Note  that  formally,  Eqs.  (8)-(10)  are  only  valid  for 
dR/dz  «  7,  so  that  results  near  the  rod  tip  are  suspect. 

Figure  6  shows  the  applied  axial  variations  of  line  den¬ 
sity,  beam  heating  and  return  current,  and  the  computed 
R(z)  at  110  ns  from  the  solution  of  Eqs.  (9)  and  (10).  For 
z  <  0.6  cm,  beam  heating  of  the  rod  plasma  dominates 
over  ohmic  heating,  and  magnetic  pressure  is  much  less 
than  kinetic  pressure,  so  that  the  zero-dimensional  model 
can  be  employed.  For  larger  values  of  z,  axial-return 
current  has  accumulated  to  high  values  and  temperature  is 
lower  due  to  reduced  beam  deposition  and  higher  mass. 
Under  these  conditions,  ohmic  heating  dominates  and 
magnetic  pressure  retards  plasma  expansion,  so  that  the 
return-current  effects  in  the  one-dimensional  model  are 
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Figure  6.  One-dimensional  model-input  variations  vs 
axial  position  and  resulting  radius  at  1 10  ns. 

important.  These  effects  result  in  the  observed  sharp 
Schlieren-boundary  transition  between  the  two  regions. 

Figure  7  shows  a  shadow  of  the  initial  rod  superim¬ 
posed  on  Schlieren  and  interferometry  images  at  110  ns, 
the  end  of  the  power  pulse.  [6,7]  In  order  to  relate  R(z) 
from  the  one-dimensional  model  (Fig.  6)  to  the  measure¬ 
ments,  the  electron  density  at  the  Schlieren  boundary  has 
been  determined  by  fringe  counting  and  Abel  inversion. 
This  electron  density,  about  4xl018  cm"3,  is  set  equal  to 
Zn(r,R)  using  computed  Z  and  R  values  at  110  ns  and 
various  axial  positions.  Solving  for  r  then  yields  the 
computed  Schlieren  boundary  as  a  function  of  z.  The 
results  of  this  calculation  are  also  shown  in  Fig.  7.  Fair 
agreement  between  the  measured  and  computed  Schlieren 
boundaries  demonstrate  the  importance  of  return-current 
and  magnetic-field  effects  to  the  rod-plasma  expansion 
dynamics  behind  the  rod  tip. 


Figure  7.  Schlieren,  interferometry,  initial-rod  images, 
and  computed  Schlieren  boundary  (circles)  at 
110  ns. 


In  conclusion,  measured  Schlieren  images  and  line- 
spread  distributions  compare  well  with  simple  zero-  and 
one-dimensional  self-similar  MHD  modeling  of  tungsten- 
plasma  expansion  for  the  plasma-filled  rod  pinch.  Close 
to  the  tip,  beam  heating  of  the  rod  plasma  and  kinetic 
pressure  dominate  the  dynamics  so  that  the  zero-dimen¬ 
sional  model  can  be  employed  to  estimate  how  expansion 
alters  the  radiographic  source.  Behind  the  tip,  accumula¬ 
ted  axial-retum  current  and  reduced  beam  deposition  lead 
to  high  ohmic  heating  and  a  magnetic  pressure  sufficient 
to  retard  plasma  expansion,  so  that  the  more  sophisticated 
analysis  is  required. 

Future  plans  include  2-dimensional  MHD  and  particle- 
in-cell  (PIC)  simulations  of  the  PFRP.  Challenges  for  this 
objective  will  be  to  model:  rod  retum-current-heating 
effects  during  the  run-down  phase,  effects  of  adsorbed 
gases  in  the  rod,  a  self-consistent  transition  between  run¬ 
down  and  plasma-opening  phases,  and  beam-  and  plasma- 
current  distributions  in  the  expanding-rod  plasma  after 
opening.  It  is  hoped  that  additional  modeling  will  help  in 
the  development  of  geometries  that  reduce  the  MHD 
wings  in  the  line  spread  to  enhance  the  radiographic  appli¬ 
cations  of  the  plasma-filled  rod  pinch. 
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